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Abstract 

We provide general determinant formulae for all n-particle form factors 
related to the trace of the energy momentum tensor and the analogue of the 
order and disorder operator in the S'f7(3)2-homogeneous Sine-Gordon model. 
We employ the form factors related to the trace of the energy momentum 
tensor in the application of the c-theorem and find perfect agreement with 
the physical picture recently obtained by means of the thermodynamic Bethe 
ansatz. For finite resonance parameter we recover the expected WZNW-coset 
central charge and for infinite resonance parameter the theory decouples into 
two free fermions. 
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1 Introduction 



In |T[] a certain physical picture for the quantum field theory of the Homogeneous 
Sine-Gordon models (HSG) was extracted from a thermodynamic Bethe ansatz 
analysis. The central aim of this manuscript is to inspect the picture for consistency 
by means of the form factor approach ^. 

The HSG-models have been constructed as integrable perturbations of WZNW- 
models. The related scattering matrices belong to a general class H, which 
describe the scattering of particles labeled by two quantum numbers, where each of 
them may be associated to a simple Lie algebra. Characteristic features of these S- 
matrices are the breaking of the parity invariance of some amplitudes and in addition 
the presence of a resonance parameter which enables the formation of unstable 
bound states. In |I| we recovered the expected Virasoro coset central charge and 
found that when the resonance parameter tends to infinity the system decouples 
into several copies of minimal affine Toda field theories. Since the ultraviolet central 
charge is also accessible by the c-theorem, the findings in may be checked for 
consistency. 

Our manuscript is organised as follows: In section 2 we recall the general prop- 
erties of form factors. In section 3 we specialise the equations to the case of the 
S't/(3)2-HSG model and provide the general solutions related to the energy momen- 
tum tensor and the analogue of the order and disorder operators. Our conclusions 
and a further outlook are presented in section 4. 

2 Generalities on Form Factors 

In order to fix our conventions and to set up the general framework we commence by 
recalling briefiy some general properties of form factors. For a proper justification of 
them in terms of general principles of quantum field theory and analytic properties 
in the complex plane we refer the reader to P, ^ 

Form factors are tensor valued functions, representing matrix elements of some 
local operator 0{x) at the origin between a multiparticle in-state and the vacuum, 
which we denote by 

F^\^^-^-{e„ ...,e^):= {o\omv,^{e,)v,^id,) . . . i^,„(^„))^^ . (i) 

Here the V^(6') are some vertex operators representing a particle of species /i de- 
pending on the rapidity 6 satisfying the so-called Zamolodchikov algebra. 

As a consequence of CPT-invariance or the braiding of two operators V^{6) one 
obtains 

T-iC'l---A',Mi+i---/ a a \ _ pC'l-'-Mi+iMi--- / n n \c ( \ (0\ 

-Tn [. . . . . .) — Pn [. . . ,Uij^i,Ui, . . .)afj,^f_,.^^[Oi^i+i) . [Z) 
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As usual we abbreviate 6ij = Oi — 6j. The analytic continuation in the complex 
^-plane at the cuts when 6 = 2Tii together with crossing leads to 

+ 2vr«, . . . , = l'^-^"'^i(^2, . . . , ^n, ^i) . (3) 

Since we are describing relativistically invariant theories we expect for an operator 
O with spin s 

+ A, . . . , ^„ + A) = e^^Ff I'^i-'^" 



'1, 



, Or 



(4) 



For a form factor whose first two particles are conjugate to each other we have a 
kinematical pole at ivr, which leads to a recursive equation relating the (n-2)- and 
the n-particle form factor 



_Res Ffif^-^" 



'1' 



1=1 



5 On)^ 



(5) 



with u being the factor of local commutativity and Jx the anti-particle of /i. We 
restrict our initial considerations to a model in which stable bound states may not 
be formed and therefore we do not need to report the so-called bound state residue 
equation. 

To be able to associate a solution of the equations (|)-(|) to a particular oper- 
ator, the following upper bound on the asymptotic behaviour [0] 



pO\^l^...|J.^l 



'1, 



, Or. 



< A 



(6) 



turns out to be very useful. Here A denotes the conformal dimension of the operator 
O in the conformal limit. For convenience we introduced the short hand notation 



limfi 



f{0. 



.On) = const exp([/(6' 



1) 



, On)]iOi) 



Ultimately form factors serve to compute correlation functions, but they may 
also be exploited to extract various other properties as for instance the difference 
between the ultraviolet and infrared Virasoro central charges, as stated in the so- 
called c-theorem |Ty 



9 



n=l Ml---Mr 



n\{2 



71 



dOi . . . dOn 



ELi"^/., cosh^, 



pO\/ii...fi„i 



On) ■ (7) 



It is essentially the property (|^) which we wish to employ for our purposes and 
check for consistency of the physical picture which emerged in [|l|. 



3 The SU(3)2-HSG model 

For finite resonance parameter the SU(3)2-HSG model describes the WZNW-coset 
model with central charge c = 6/5 perturbed by an operator with conformal di- 
mension A = 3/5. The model contains only two self-conjugate solitons which are 
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conveniently denoted by "+" and "— ", since that will allow for compact notations. 
The S-matrix elements read 

1 / TT 

S±± = —1 and S-i-^{9) = ±tanh - [9 ± a — i— 

This means the scattering of particles of the same type is simply described by the 
S-matrix of the thermal perturbation of the Ising model. Also the remaining am- 
plitudes do not possess poles inside the physical sheet, such that the formation of 
stable particles via fusing is not possible. For vanishing resonance parameter a the 
amplitudes 5"^^= coincides formally with the ones which describe the massless flow 



between the tricritical Ising and the critical Ising model as analysed in []TT|. How- 
ever, there is an important conceptual difference since we view the expressions (H) 
as describing the scattering of massive particles. This has important consequences 
on the construction of the form factors and in fact the solution we compute below 
will be different from the one proposed in []TT|. In the HSG setting the massless flow 
was recovered in the context of the thermodynamic Bethe ansatz [|l| only as a sub- 
system in terms of specially introduced variables combining the inverse temperature 
and the resonance parameter. When the resonance parameter tends to infinity the 
amplitudes S±^ become one, describing non-interacting scattering, such that the 
-system and the "-"-system decouple. 

Attempting now to solve the equations presented in section 2, we proceed as 
usual in this context and we make a factorization ansatz which already extracts 
explicitly some of the singularity structure we expect to find. For the case at hand 
we have to have a kinematical pole at in when two particles are conjugate to each 
other 



I X ± m X =F 



F:x'i9^,) 



(9) 

We introduced the variable Xi = exp9i. The are normalization constants. 

As common we suppose that the so-called minimal form factor satisfies 

Fl^iO) = ^n(-^^)^.(^) = i^n(2vr^ - 9) (10) 

and has neither zeros nor poles in the physical sheet. Then, if we further assume that 
Qn^^^"'^"{9i, . . . , 9n) is separately symmetric in the first / and the last m rapidities 
and in addition 27r2-periodic function in all rapidities, the ansatz solves Watson's 
equations (H) and @ by construction. In particular we have 

I X + m X — mx— ;x4- 



Vn \Xi, . . . ,Xn) — Vn \Xl+l, ■ ■ ■ X„, Xi, . . . X/J 



'111 
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such that when we have constructed a solution for one particular ordering of the 
/i's, e.g. the upper sign in @, we can obtain the solution for a permuted ordering 
by the monodromy properties. Especially the reversed order we obtain by apply- 
ing equation (|lT]). Despite the fact that we do not gain anything new, it is still 
instructive to verify as a consistency check also for the different ordering. The 
monodromy properties allow some simplification in the notation and from now on 
we restrict our attention w.l.g. to the upper sign in (P). In addition we deduce 
from equation that for a spinless operator O the total degree of Q'^ has to be 
/(/ - l)/2- m(m - l)/2. 

A solution for the minimal form factors, i.e. of equations is found easily 



—i sinh - 




A/'^(^)exp 



(12) 
(13) 

(14) 



Here F^:^^{6) is the well-known minimal form factor of the thermally perturbed 
Ising model |T! 

for vanishing a up to normalization with the expression found in |]lT 
the normalization function M^{9) = exp |^!Z!li±iM^ 



1^ and for the upper choice of the signs, equation (|T^) coincides 

^ ^ We introduced 



, with G = 0.91597 being 

the Catalan constant. The minimal form factors possess various properties which we 
would like to employ in the course of our argumentation. They obey the functional 
identities 



mm V 



— sinh 6 
2 

. 2=F1 



4-S. 

12 6 2 



cosh \{e±a-'^ 
We will also exploit the asymptotic behaviour 

limF±i^(±^)~e-^, 



(15) 
(16) 



1 

2' 




(17) 



Together with the factorization ansatz (P) this leads us immediately to the relations 





i 




+ 

i 




i 




+ 

i 



m — I — 1 



for l<i<l 
for I < i < n . 



(18) 
(19) 
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which are useful in the identification process of a particular solution with a specific 
operator. Since we may restrict our attention to one particular ordering only, we 
abbreviate the r.h.s. of (|^) from now on as F^'''*" and similar for the Q's. 

Substituting the ansatz (P) into the kinematic residue equation (|^) reduces, with 
the help of (p!5| ) and (|T6|) , the whole problem of determining the form factors to the 
following recursive equations 

Qn+2 X, . . . , Xn) = D^i (^1 ; • • • ; ^n)Qn ' ' {^lj---j^n) (20) 

-| m 

= -Hxy+V+EHe-x)-'=(l-c.(-iy+V. (21) 

^ k=0 

Here we introduced yet another short hand notation, namely for elementary sym- 
metric polynomials a^ixi, . . . ,xi) = and (Tfc(x/+i, . . . , x„) = cr^ @. Below we 
shall also employ ak when the polynomials depend on all n variables, ak when 
they depend on the n inverse variables, i.e. a;~^and ak when they depend on the n 
variables Xie~'^. 

The recursive equations for the constants turn out to be 

Fixing one of the lowest constants, the solutions to these equations read 

j^O\2s+t,m j^s'm2^s{2s—m~l+2t)^s'ma/2j^O\t,m ^ 1 (23) 

For specific operators we will provide below the explicit expressions for the ff'^l''™-. 
Notice that there is a certain ambiguity contained in the equations (|2^), i.e. we can 
multiply if^l^'"* by i^^, P'' or (—1)' and produce a new solution. However, since in 
practical applications we are usually dealing with the absolute values of the form 
factors, these ambiguities will turn out to be irrelevant. 



3.1 Solutions 

Whenever we consider F^'^'™ with / even for vanishing resonance parameter a, we 
can use the kinematic residue equation //2-times and finally construct F^'"''", 
which should correspond to a form factor of the thermally perturbed Ising model. 
In other words in that case we can always use the well-known solutions as 
the initial condition for the recursive problem (pO|). 

*The elementary symmetric polynomials are generated by 



n n 

+ Xk) = ^a;""''crfe(xi, . . . ,x„) , i.e. 

fc=l k=0 

(For more properties see e.g. iH.) 



1 f dz 
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3.1.1 The energy momentum tensor O 



The only non- vanishing form factor of the energy momentum tensor in the thermaUy 
perturbed Ising model is well know to be 



F®{e) = -2mm^smh{e/2) 



(24) 



From this equation we deduce immediately that [-F®l''^]j = 1/2, which serves on the 
other hand to fix with the help of (|1|) and ([T|). Recalling that the energy 

momentum tensor is proportional to the perturbing field and the fact that the 
conformal dimension of the latter is A = 3/5 for the 5'f/(2)3-HSG model, the value 

= 1/2 is compatible with the bound (||). As a further consequence of 
we deduce 



H 



e|o,2 



27rm^ 



(25) 



as the initial value for the computation of all higher constants in ([23|). The distinc- 
tion between m_ and m+ indicates that in principle the mass scales could be very 
different as discussed in Notice that is reached only formally, since the 

kinematic residue equation does not connect to the vacuum expectation value. The 
initial values for the recursive equations (EH) are 



^e|o,2 

V2 



and 



^e|o,2t 

V2t 



for t > 2 . 



(26) 



Taking now uj = 1, the solutions to (|20|) , with the same asymptotic behaviour as the 
energy momentum tensor in the thermally perturbed Ising model, are computed to 



Qf':i' = (-l)(^+^)*e-*Viai(a+)^-*(a2-J^-Met^« for t > 1, s > 1, 



where A® is a (t + s — 2)x (t + s — 2)-matrix whose entries are given by 



2(j-i)+l 



Explicitly we have 



^3 



-^5 

(In 



On 



^7 



for \<i<t 
for t < 2 < s + t 



\ 




(27) 



f28) 








(29) 



(-l)*'^2Vi / 
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One may easily verify case-by-case that ( p?] ) is a solution of (^2j) to relatively high 



orders in s and t. A general proof of this result, which we present elsewhere ||16|| , 
can be obtained by exploiting the fact that the determinant of A may also be 
represented in terms contour integrals 

2s t— 1 

det^® = j^du,...f dut^^ fdv,...f dvs-i n n (30) 

1=1 j=l Uj 

2s+2t s-l „, I - s-lt-1 

X n n;%iA n n (-i--nnn(-?+-j)- 

i=l+2sj=l l<«<i<t-l l<i<j<s-l j=li=l 

In order to establish the equivalence between (^) and ( pO] ) we simply use the 
integral representation for the symmetric polynomals as stated in the footnote. 
The integrals in (jS^) are understood a.s § dz = {2TTi)^^ f^^^^^dz with g being an 
arbitrary positive real number. 

Assembling now all the quantities we obtain for instance 

P 1 - -^^-e^'""^'"^^'(2 + E.<,cosh(%)) 

F, (^1, ^„ ^3, ^4) - 2cOsh(ei,/2)cOsh(034/2) Q ^n.in (^.0 • (31) 

Having computed all form factors for the energy momentum tensor we are in the 
position to apply the c-theorem, i.e. we can in principle evaluate (|^. For finite 
values of a we obtain 

Ac(2) = 1, Ac(^) = 1.197..., Ac(^) = 1.199 . . . , for a < oo (32) 

where in the notation Ac^^\ the superscript n indicates the upper limit in (^. 
Thus, the expected value of c = 6/5 is well reproduced. Apart from Ac*^^-*, in which 
case the calculation can be performed anal5d;ically, the integrals in (0) are computed 
directly by a brute force Monte Carlo integration. 

When the resonance parameter tends to infinity the system decouples and we 
are left with two non-interacting free fermions, such that the only contribution in 
the sum is twice the free fermion two-particle contribution, such that 

lim Ac = 1 . (33) 
In order to see this we collect the leading order behaviours form our general solution 



Jim i/SSf ~ e-^ Jim Q^T ~ e'^'^^-'^% Jim H F^;^ (6,) ~ e--^ (34) 

i<j 



which means 

lim fS^' ~ e-(*+^-i)'^ . (35) 



tr— >oo 



Hence the only non- vanishing form factors in this limit are F2 ' and F2 ' , which 
establishes (^3]). 
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3.1.2 The order operator S 

For the other sectors we may proceed similarly, i.e. viewing always the thermally 
perturbed Ising model as a benchmark. Taking now u = 1, we recall the solution 
for the order operator 



-^2^ + 1(^15 • • • 5 ^2s + ly 



Kj 



(36) 



With this information we may fix the initial values of the recursive equations (|20|) 
and (HTi at once to 



Q 



2t+l 



ia,,^,Y and H^^^'^ = Ff . 



(37) 



Furthermore, we deduce from equation (^) that = 0. Respecting these 

constraints we find as explicit solutions 



Q2i?2m = (-l)^^+'^*(^i)^(4J^-*-'(^r)-^Km)~*det^^ 
where J\P is a (t + s) x (t + s)-matrix whose entries are given by 



Explicitly this reads 
/ 



a. 



2(i-i) 



V 



(j-i- 






-1 


1 


at 









1 


at 
















^1 













^3 

















for l<i<t 
for t < i < s + t 



\ 




0"2s 




(-l)(*+^)^2m ) 



(3J 



(39) 



(40) 



Once again the determinant of A admits an integral representation 



2s t 



detA^ = i-iy'^'-^^ j^dui...j^dut j^dvi...j^dv, fill 



2s+2t+l s „, I ^. 

X n niSA n 



Uj J I fXj 2 

i=l j=l Uj 

s t 



(41) 



i=l+2s j=l 



/ 2 2 



l<i<j<i 



' n (-j--f)nn(-?+-i) 

l<i<jr<s j=l i=l 



which may be used for a general proof [1^ . 
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When the resonance parameter tends to infinity we obtain the following asymp- 
totic behaviour 

jim Q£^2m ~ e-'^ (42) 
;imi7£S'n^mm^(%) = ^onst H F^^{9,j) [] ^mi^l^) • (43) 

i<j l<i<j<2s 2s<i<j<2s+2t+l 

This means unless s = 0, that is a reduction to the thermally perturbed Ising model, 
the form factors will vanish in this limit. 



3.1.3 The disorder operator n 



For the disorder operator we have o; = — 1 and the solution acquires the same form 
as in the previous case 



i<j 



(44) 



Similar as for the order variable we can fix the initial values of the recursive equa- 



tions dg) and (H) to 

Q2i ' = [cr2t) 

Furthermore, we deduce [F^'^"*'^*] 
general solution 

\sti 



[a2t) ' and = F^^ . (45) 

0. Respecting these constraints we find 



Qfef = {-ir{a,,^2t)^-\alr'ia^r'detA^. 
where is a {t + s)x[t + s)-matrix whose entries are given by 



Explicitly we have 



A'' 



'2(i-i) 
.(-l)(^--+*)a-(^._,+,) 



/ 1 




1 





(^2 


(^4 




1 


at 













^2 






1 


-^2 














for l<i<t 
for t < i < s + t 



\ 




(^2s 






(46) 



(47) 



(4J 



■^y^2t J 
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Similarly as in the previous sections we can write the determinant of A alternatively 
in form of an integral representation 

2s t 

detA^ = (-l)^(^-^ idu^..Jdutidv^..Jdvsl[t[ '^'.X"^; (49) 

i=i j=i 

2s+2t s _i_ - St 

X n n n (-i--Dnn(-^+-i)- 

i=l+2sj=l'^j l<i<j<t l<i<j<s j=li=l 

When the resonance parameter tends to infinity we observe the following asymp- 
totic behaviour 

jimQ^^^f = (-l)^*gr'°Qf''* (50) 
Jim/ZriS* n^mm^(%) = const n F++(%) [] ^mi;(%) (51) 

i<j l<i<j<2s 2s<i<j<2t+2s 

such that 

\imF,f:i'-F,f'''F,T'r (52) 

This means also in this sector we observe the decoupling of the theory into two free 
fermions. 



4 Conclusions 

The application of the c-theorem confirms very well the physical picture we found 
in |]T[] from the thermodynamic Bethe ansatz. For finite resonance parameter we 
recover the expected Virasoro central charge of c = 6/5 and for a ^ oo the theory 
decouples in all sectors into two non-interacting free fermions. Besides the con- 
struction of all n-particle form factors related to the trace of energy momentum, 
we computed in addition the complete solutions for the order and disorder operator 
in form of determinants whose entries are symmetric polynomials. Such determi- 
nant formulae have occurred before in various places in the literature, e.g. |]^, |T7 



Representing the solutions for form factors in this form has turned out to be use- 



ful in the construction of correlation functions and might eventually lead to a 
reformulation of the whole problem in terms of differential equations analogous to 
the situation in conformal field theory |jl9| . Apart from higher spin solutions which 



may always be constructed by including the polynomials as suggested in ^0[, we 
did not find any additional solutions related to other sectors. We expect that a 
careful analysis of the cluster decomposition property will lead to more conclusive 
statements concerning the question whether such solutions exist at all. From a 
mathematical point of view it is also desirable to present a rigorous proof of the 



determinant formulae \16 
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